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Abstract 

Using first principles, it is demonstrated that radiative damping alone cannot 
lead to a nonvanishing electro-optic effect in a chiral isotropic medium. This 
conclusion is in contrast with that obtained by a calculation in which damping 
effects are included using the standard phenomenological model. We show 
that these predictions differ because the phenomenological damping equations 
are valid only in regions where the frequencies of the applied electromagnetic 
fields are nearly resonant with the atomic transitions. We also show that 
collisional damping can lead to a nonvanishing electrooptic effect, but with a 
strength sufficiently weak that it is unlikely to be observable under realistic 
laboratory conditions. 
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Several recent papers [1,2] have discussed the question of properly taking into account 
various relaxation processes while calculating the nonlinear response of an optical system. 
Even the existence of certain nonlinear optical processes is thought to be closely linked to 
the existence of a damping mechanism [3-5]. In this connection, it is especially important 
to incorporate in a consistent manner the effects of relaxation processes. Very often the 
nonlinear response [6] is calculated by modifying the equation for the off-diagonal elements 
of the density matrix (coherences) by introducing phenomenological relaxation terms as 
follows: 

dpij 

—-— — —iuJijPij + field terms =^ (1) 



dp 
~dt 



—icUijPij — TijPij + different field terms. 



(2) 



The equations for the populations are also modified appropriately. Such modifications have 
been extensively used in nonlinear optics and even have led to the prediction of new effects 
such as collision-induced resonances which have been subsequently been observved exper- 
imentally [7]. Kauranen and Persoons [4] have recently presented a theoretical argument 
that predicts the existence of a linear electro-optic effect (EOE) in chiral isotropic media 
provided material damping is taken into account. Their result follows by using (2). However, 
it is not clear a priori if Eq. (2) can be used to describe the linear electro-optic effect. In 
order to see the origin of this uncertainty, let us examine the expression for the nonlinear 
susceptibility describing the electro-optic effect in a chiral isotropic medium. The derivation 
given in Ref. [4] is based on the standard phenomenological equations (2) which take into 
account various damping processes in the medium. The nonlinear susceptibility is shown to 
have contributions of the form 
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The authors of ref. [4] have suggested that this damping-dependent contribution is the one 
which can lead to a nonvanishing electro-optic effect in a chiral isotropic medium. Let us 
examine this contribution further. We note first that the usual expression for the second- 
order susceptibility consists of two energy denominators whereas the above contribution 
consists of three. Clearly such a term arises from the combination of two contributions as 
X can be written as 



X 



ij^mg ~l~ ^ ~l~ i^mg') J^ng "^Ing ^ng ~l~ ^l^ng _ 
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We note also that denominators such as {ung — ijng) do not have an optical frequency con- 
tribution. Such denominators arise from the interaction of the system with a zero-frequency 
field. We show below that in a correct treatment of radiative damping, the denominator 
should be replaced by ones that involve frequency-dependent damping coefficients. Thus a 
first-principles treatment would lead to 



X = 



{(Jmg + U; + nmg{i^)) 
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Note that the frequency dependence of 7 in each denominator depends on the frequency 
component of the electromagnetic field responsible for such a denominator. Thus the de- 
nominators corresponding to the static field have dampings evaluated at zero frequency. As 
discussed below, for the case of radiative damping 7ng(0) vanishes identically, which implies 
that X — 0. Thus a first principles (and correct) treatment of radiative damping does not 
lead to any electro-optic effect in a chiral isotropic medium. We also show below that X 
is at most very small for the case of collisional damping. The nonvanishing of the EOE 
effect reported earlier is due to inappropriate use of equations which are not valid for the 
calculation of the EOE effect. Thus, when using the modification (2) in the calculation of 
the nonlinear optical response, one has to keep in view the conditions under which (2) has 
been derived. This need necessitates an examination of the microscopic theory leading to the 
derivation of the result (2). It may also be noted that, in recent times, one has discovered a 
number of other interesting situations which cannot be described by equations like (2). For 
example, there are situations under which the coherences get coupled to the populations, 
and this situation has led to considerable work on quantum interferences [8]. In addition, 
there is the subject of inhibited spontaneous emission where the modifications of (2) due to 
strong external fields play an important role [9]. 

In order to uncover the role of relaxation mechanisms on the response to external fields 
and to determine how relaxation depends on the frequency of the applied field, we consider 
first a very simple model. This model brings out the salient features of the problem and 
enables us to establish that the form of the damping operator depends on the various fre- 
quency scales in the system. We consider the case in which the medium can be described 
by a one-dimensional harmonic oscillator with displacement x and with frequency Uq- Let 
the medium interact with an external electromagnetic field of frequency u described by 

E = £e-''"* + S*e'^K (6) 

The equation of motion with a phenomenological damping constant F is 

eS 

x + rx + ujIx = — e""^* + c.c. . (7) 
m 

The response of the medium can then be expressed as 

e =xM^e-^'"* + c.c., (8) 



xH = —r^ 2 — (9) 

m{U!Q — CU^ — ICUL). 

In this manner one obtains the familiar response function. We would like to examine whether 
the response x{^) given by Eq. (9) is valid for all frequencies. Thus we would like to 
understand if the introduction of a frequency-independent damping constant F in Eq. (7) is 
justified for all frequencies of the applied electromagnetic field. For this purpose we start 
from first principles. Let us consider the interaction of the system oscillator with a bath. 
The bath will be responsible for the relaxation processes described phenomenologically by 
the damping parameter F in Eq. (7). As usual we model the bath by a set of harmonic 
oscillators. The Hamiltonian for the system oscillator interacting with a bath is given by 
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H^— + -mu;lx'^-exS{t)-xF{t), (10) 

where S{t) is the time-dependent electromagnetic field and F(t) represents the effect of the 
bath terms 

n^) = E(W"^'^^* + h.c.). (11) 
j 

Here u!j{> 0) are the frequencies of the bath oscillators aj and gj are the coupling constants 
of the system oscillator with the bath oscillators. The Heisenberg equations can be easily 
derived from (10): 

x — p/m, p — —mujQX + e£{t) + F{t), 

dj ^ig*e''^^^x(t). (12) 
We integrate formally the equation for aj and substitute it into the equation for p to obtain 

p = -mouox + e£{t) + Fo(t) + /* K{t - T)x{T)dT, (13) 

Jo 

where 

^o(t) = E^j%(0)e"'"'* + h-c. , (14) 



K{t - r) = (i ^ |^,f e-^'^^(*-^) + c.c). (15) 
j 

Note that Eq. (13) is derived without any approximation. The further simplification will 
depend on the values of \gj\,ujj,uj etc. Let us examine the average response for the case in 
which S{t) = £e~^'^*' + c.c. Note that the mean value of the operator 0^(0) is zero and hence 
{Fo{t)) = 0. It should be borne in mind that u; is positive. Using Eqs. (13) and (14), taking 
quantum mechanical expectation values and the long-time limit i — > oo, we obtain 

^"^ = m(a;g-a;^)-ir(a;) +""' ^''^ 



where 



K(uj) = limT^ |(7iPi I ; 



u) e — i(oi)j + oj) 



^K\u)+iK"{u), (17) 



K"{uj)^Y.\93?^^i^3-^)- (18) 

3 

The exact result (16) has the same structure as (9) except for the important difference that 
a;r is replaced by a function K [uj) which is dependent on the frequency uj of the applied 
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electro-magnetic field. In addition there is a dispersive contribution Re K(a;). Note further 
that very often one replaces (18) by 

K"{u;)^Y.\9j\'^^i'^j-^o). (19) 
j 

Clearly this can be done if the frequency loq of the system oscillator is very close to the 
applied frequency, i.e., essentially in the resonance region. If the frequency u happens to be 
far away from a resonance frequency, then the phenomenological equation (7) should not be 
used. This is the real reason why usage of equations like (2) and (7) can give rise to incorrect 
nonlinear optical response for applied frequencies far away from the transition frequencies. 
We also find from (18) that for static response 

limi^"(u;)^7r|^,f| ^0, (20) 

for the usual radiative coupling. Thus the static response functions would be independent 
of the damping term. More generally no damping term can appear in the static response as 
long as the bath does not have a characteristic static frequency. 

The features discussed above are valid rather generally. To see this we consider the 
dynamical equations for a two-level system undergoing, say, radiative damping. The case 
of a two-level system is more involved because of the intrinsic nonlinearity of the two level 
system. However the salient features can be uncovered by using the wavefunction approach. 
Let us write the interaction Hamiltonian of a two level system interacting with the field and 
undergoing radiative damping, as 

H = ^a;o|e)(e| - h[G{t)\e){g\ + h.c] - hJ2{gkake-''^'^'\e){g\ + h.c), (21) 

k 

where we sum over all field modes labelled by the index k and where, 

G{t) = d.E{t)/h. (22) 
= Goe"^"^* + c.c. . 

The last term in (21) is responsible for the radiative decay of the atom. The couphng to the 
mode k with frequency uo^. of the electromagnetic field is represented by g^^ is the photon 
annihilation operator. The wavefunction of the whole system can be expressed as 

\^) = Ve|e, {0}) + ^Pg\g, {0}) + Y.^k\9, {k}), (23) 

k 

where {0}({A;}) represents the vacuum (one photon in mode k) state of the field. The 
Schrodinger equation leads to 

ijg = tG*{t)iJe, 

k 

i^k = igli^ee''^''. (24) 
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The initial conditions are i/jg = i^k = = 1- The induced polarization is to be ob- 

tained from the off-diagonal element Peg — '0e'0p- Note that to first order in the applied 
electromagnetic field, p^g is 

p^:]{t)^€\t)rj~'\t)+i^f\t)i^f\t) 

-€\t) (25) 

tl2^ ^(+)e-^'-* + V;(-)e''^*. (26) 
To obtain the steady state response we combine last two equations in (24) 

4 = -^a;o^e + iG{t)^g - J2 \9k\' f e-*"'=(*-^Ve(T)ciT, (27) 

k 

and thus to first order in the external electromagnetic field we obtain 

i^f^ = -iu;oV'l'^ + iG{t) - Y: \gk\' f e--'=(*--)V'^^)(T)dT. (28) 

k 

In terms of Laplace transforms we have the result 

where we have used the explicit form (22) From (29) we get the response in the long time 
limit 

-0^+) = -uj) + Y^ \9k?ii^k - i^)"^ Go, 

= (^i{uo+u;) + ^\gk\\iUk + iu)~'^ G*. (30) 
The induced polarization can now be calculated 

p{t) = {Pegdlg + C.C.) 

= Poe-'^' + c.c. , (31) 

where po is calculated using Eq. (30) as 

Po = V^+)4 + Vi-^*4. (32) 

This is the most general result for the linear response. No assumption has been made 
regarding the nature of the bath. It should be borne in mind that all frequencies in (30) are 
positive. The radiative corrections enter the response function through the quantity 

K{z)=J2\gk\\z + zu;k)-'. (33) 

k 
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It should be noted that the actual radiative correction terms depend on the frequencies of 
the applied fields rather than the atomic frequencies. It is only when the applied frequency 
is close to the atomic frequency that we can use the approximate replacement a; — > a;o in 
■0^+) (this cannot be done in ijj^''^). We thus find that the counter- rotating contribution ■0^'")* 
in Eq. (32) does not depend on the radiative damping [10]. The rotating- wave contribution 
depends on the radiative damping; however the radiative damping is to be cvahiatcd at the 
applied frequency. If such an applied frequency is very far from the atomic transition, as for 
example for dc fields, then no radiative damping term appears in the response. Thus the 
full quantum mechanical calculation also leads to the same conclusion as we derived for the 
simple oscillator model. Further the above analysis can be easily extended to the multilevel 
systems and to the calculation of second-order and higher-order response. We find similar 
conclusions regarding the various denominators which appear in response functions. The 
argument given in the context of Eqs. (4) and (5) is correct and we rule out the possibility 
of the occurence of electro-optics effect due to radiative damping. 

A pertinent question could be: can other damping mechanisms, such as phase changing 
collisions, possibly lead to the nonvanishing of the EOE in isotropic chiral medium? This 
question has to be examined by considering a detailed microscopic model for the collisional 
process. However, a simple model calculation outlined below suggests that even if the effect 
is non-vanishing it must be extremely small; in particular it must exponentially small in a 
large quantity. 

Consider the equation for the optical coherence a = peg- Let /(t) be a stochestic source 
which represents the effect of phase changing collisions. We model f{t) to be a Gaussian 
stochastic process with correlations given by 

(fit)) = 0, if it) fir)) = e-^l*-^l fl (34) 

Here F"^ is the magnitude of the collision time. The equation for the optical coherence can 
be written in the form 



&^-iAa-if{t)a + iG, (35) 

where G represents the external field. If F~^ is the smallest time scale in the problem, then 
one can show using the standard methods [11] that 

M= (f + ^Aj \iG). (36) 

Thus one recovers the result of the phenomenological theory. However, for the response to 
a static field, A is of the order of the optical frequency whereas typical collisional process 
take place over a scale which is of the order of psec or larger. Thus F^^ is no longer the 
smallest time in the problem. The smallest time scale will instead be A~^. In such a case 
one can show that in the long-time limit 




^iG r dTe-'^^e--^fo^\ (37) 
Jo 
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Note that the square bracket in Eq. (4) is just the real part of dre '"'"s'^, and thus 

if we had treated the damping properly it has to be replaced by 



Thus coUisional damping can make the EOE in chiral isotropic medium nonzero. However it 
would be extremely small unless the strength of collisions is comparable to ujny i.e., ujng ~ fng- 

In conclusion, we have shown that radiative damping cannot lead to a nonvanishing EOE 
in a chiral isotropic material. For the case of coUisional damping, a nonvanishing EOE is 
predicted, but the magnitude of this effect is expected to be so small that it is unlikely 
that this effect could be observed experimentally. These results are in contrast with recent 
suggestions that relaxation effects can lead to an EOE in chiral isotropic materials, with 
potentially important practical implications. More generally, we have shown that in general 
it is not adequate to use a frequency-independent damping parameter in treating relaxation 
processes within the context of density matrix calculations. 
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